Renormalization Group Theory of the Three-Dimensional Dilute 

Bose Gas 



M. Bijlsma and H.T.C. Stoof 
Institute of Theoretical Physics, University of Utrecht, Princetonplein 5 
P.O. Box 80.006, 3508 TA Utrecht, The Netherlands 

Abstract 

We study the three-dimensional atomic Bose gas using renormahzation group 
techniques. Using our knowledge of the microscopic details of the interatomic 
interaction, we determine the correct initial values of our renormalization 
group equations and thus obtain also information on nonuniversal properties. 
As a result, we can predict for instance the critical temperature of the gas and 
the superfluid and condensate density of the Bose-Einstein condensed phase 
in the regime naAj^ <C 1. 
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I. INTRODUCTION 



After a long history in which a large number of experimental groups around the world 
contributed to the development of succesful methods to master stabilization and cooling 
of dilute Bose gases, last year the aim of achieving Bose-Einstein condensation in such a 
system was finally reached. Indeed, a macroscopic occupation of the one-particle ground 
state was irrefutably observed in magnetically trapped and evaporatively cooled alkali gas 
samples of ^''Rb and ^^Na using relatively simple time-of-fiight measurements The 
transition that was claimed to be seen in an experiment using '^Li was less convincing [Q. 
In the latter case, the interatomic interaction is effectively attractive and the potential 
has a negative scattering length a. Therefore, Bose-Einstein condensation in this system is 
preempted by a first order phase transition to a liquid or solid phase in the homogeneous case 
Nevertheless, for inhomogeneous gas samples the trapping potential has a stabilizing 
influence and a macroscopic occupation of the ground state is possible in principle. However, 
when the condensate contains more than some 1500 particles under the conditions of the 
''Li experiment 0, the condensate is still expected to collapse Q. 

After these first experiments, which were primarily aimed at proving the existence of a 
Bose condensate, many experimental groups are now building or improving on their experi- 
mental setups to be able to perform much more precise measurements of various interesting 
properties of the gas in the degenerate regime. Superfluidity the condensate density 
and its profile, the dynamics of condensate formation p|, the Josephson effect [§, vortex 
dynamics, collective excitations [[I0l -[l3[| and the precise value of the critical temperature are 
examples of phenomena and quantities of interest. Other types of experiments will presum- 
ably also study the properties of mixtures of atomic gases. In this respect one might think of 
two bosonic species with a different sign of the scattering length, or mixtures of bosons and 
fermions, or mixtures uniting both aspects. In the case of a pure fermionic gas of ^Li atoms, 
a BCS transition to a superfluid state is predicted to occur and should be within reach of 
the current experimental technology |T^. Furthermore, in some cases the magnitude and 
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even the sign of the scattering length can be changed by varying the apphed magnetic field. 
This opens the road to yet another type of experiment. 

It is clear from these possibilities that a large number of experiments are expected to be 
performed in the near future, which clearly makes the degenerate dilute Bose gas also a very 
interesting subject for theoretical studies. Indeed, this field of research has rapidly expanded 
during the last year. However, most approaches to the dilute Bose gas use the Bogoliubov (or 
Popov) theory and are therefore of mean-field type and susceptible to improvements, both 
from a practical as well as a fundamental point of view. In these approaches one mostly uses 
the so-called two-body T— matrix, or equivalently the scattering length a. Technically, this 
important quantity describes the collisions taking place in the dilute Bose gas by summing all 
possible two-body scattering processes, i.e. all ladder diagrams, without taking into account 
the fact that the surrounding gaseous medium has an effect on these collisions. However, we 
have recently shown that the many-body corrections arising from the surrounding gas may 
be important ||15|, and are even essential for solving the problem connected to the order of 
the phase transition which is found to be of first order when using the two-body T— matrix 

Including quantatively the same many-body corrections in the case of a highly inhomo- 
geneous gas sample has at this point not yet been done. Moreover, introducing the effect of 
the medium on two-particle collisions also in the condensed phase by means of the many- 
body T— matrix leads to fundamental problems if we want to describe the physics at long 
wavelengths correctly as the ladder diagrams contain infrared divergencies in this case. Us- 
ing renormalization group techniques, we expect in principle to be capable of resolving these 
infrared problems as with this method a correct resummation of diagrams is automatically 
performed, eliminating any potentially troublesome large distance behavior of the individ- 
ual diagrams. Furthermore, a renormalization group calculation can be used to improve 
the usual mean-field approaches and the many-body T— matrix theory in the critical region. 
Indeed, we recently predicted by these means for example that the critical temperature in 
the ^^Rb and ^^Na experiments can, due to interaction effects only, be raised with as much 
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as 10 % compared to the ideal gas value found from the criterion Ucl^^h = C(3/2) ^ 2.612 
PCI|| , which is identical to the criterion also found using mean-field calculations |T^. Here n 
is the density and Is^th = {2Txfi^ /mkBTy/'^ the thermal de Broglie wavelength of the atoms in 
the gas. Therefore, it is very well conceivable that in the critical region also other properties 
of the dilute Bose gas, such as the superfluid and condensate densities, will significantly 
change when going beyond the mean-field level or beyond the many-body T— matrix theory. 

The renormalization group method is a very powerful method which was in first instance 
developed by Wilson ^l[ to study the universal properties of second order phase transitions. 



The basic idea is to perform the trace in the grand canonical partition function Zgr = 
j^^(^Q-P{H-'fj.j\f)^^ gradually, starting with the high momentum states. After each step one 
tries to find a new effective Hamiltonian such that Zg,. = Tr'{e^^^'^'^^'^'>) and the trace is 
limited to the low momentum states which have not been reached yet. One proceeds until 
the complete sum has been performed. Besides the partition function one in this manner 
also ends up with the effective Hamiltonian describing the long distance properties of the 
system. 

The renormalization group method has been applied to the dilute Bose gas |2^-^| before. 



but without actually performing an extensive quantitative study of this system. This is 
due to the fact that in general the quantities of nonuniversal nature, such as the critical 
temperature and the superfluid and condensate densities, depend on the microscopic details 
of the system considered. Put differently, they depend on the ultraviolet cutoff A of the 
theory and this quantity is usually unknown. However, due to the diluteness of the gas the 
nonuniversal properties are in the present and forthcoming experiments the most interesting 
ones, and therefore we are in this paper mainly interested in these aspects. The reason 
that it is nevertheless possible to perform a quantitative study of the dilute Bose gas using 
the renormalization group method, is that for this system we do have sufficient information 
about the microscopic details to calculate and predict the nonuniversal properties. Thus 
we can, by correctly applying this knowledge, use the renormalization group method and 
eliminate the cutoff dependence at the same time. We will come back to this point in Sec. 



II. Furthermore, we will show that, in contrast to regular perturbation theory, the problems 
related to the infrared divergencies are in principle indeed resolved, but lead nevertheless 
to some problems whose solution requires further investigation. However, these problems 
are only of importance when the interaction energy is no longer negligible compared to the 
kinetic energy of the particles. The dimensionless parameter reflecting this aspect is naA^^. 
Therefore, we will in this paper first concentrate on the regime where naK^^ <C 1. 

We treat here only the homogeneous Bose gas with effectively repulsive interactions, i.e. 
with a positive scattering length. However, as in all experiments up till now the number 
of particles N is so large that the critical temperature Tc is much larger than the energy 
splitting fiuj between subsequent levels, one can practically for all temperatures use a local 
density approximation to describe the gas in the trap. The criterion for this description 
to be valid is that the correlation length should be smaller than the typical length scale on 
which the atomic density varies. Therefore, a local density approximation breaks down close 
to the edge of the gas cloud, which is for most practical purposes an unimportant region, 
but also in the center of the trap if the temperature approaches the critical temperature and 
the diverging correlation length starts to exceed the typical dimensions of the trap. The 
temperature interval where this occurs has a width of 0{Tc{fiuj / kBTc)) around the critical 
temperature. As this region is very small, we conclude that the results we find in this paper 
for the homogeneous gas are essentially also valid for the inner part of the trapped Bose 
gases, and in particular pertain to the ^^Rb and ^^Na experiments. Furthermore, we want 
to remark here that it would in principle also be possible to set up a renormalization group 
calculation for the inhomogeneous Bose gas. Of course, there is no real second order phase 
transition present in this system because the correlation length ^ can never become infinite, 
but the techniques of renormalization group as presented in this article can still be used to 
calculate the (nonuniversal) properties also in this case. 

The paper is organized as follows. In Sec. II we briefly discuss the theoretical description 
of the dilute Bose gas and the renormalization group method. In Sec. HI we first apply the 
renormalization group to the uncondensed phase of the Bose gas because the flow equations 
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are relatively simple and easy to understand in this case. In Sec. IV we then go over to the 
degenerate Bose gas and again describe the gas properties following from the renormalization 
group approach. Finally, in Sec. V we end with some concluding remarks. In the numerical 
calculations we always use ^^Na as an example, because the experiment with this atomic 
species is closest to the conditions of homogeneity 0. We take in these cases the most 
up-to-date value of 52ao for the two-body scattering length p6 |. 



II. THE RENORMALIZATION GROUP 



The renormalization group equations are most easily obtained using the functional inte- 



gral formulation of the grand canonical partition function |2^. We thus write 



Z,, = Tr(e-^(«-'^)) = J dirW] exp . (1) 

This functional integral is over c-number fields '?/'*(x, r) and '?/'(x, r) periodic in imaginary 
time over h/3 = h/ksT. The so-called Euclidian action for the dilute Bose gas is given by 



h— yU 

OT 2m 



V'(x,r) 



+ ^ I f/x| rfx' V^*(x, T)ri^', r)V(x - x')V^(x', r)^(x, r)) , (2) 

with /i the chemical potential and V^(x — x') the effectively repulsive interaction potential. 
In principle, the action also contains a term describing three-particle interactions, i.e. 



and terms describing interactions between four and more particles due to the finite extent 
of the electron clouds of the atoms. However, as we are describing the dilute Bose gas these 
terms are expected to make in general no significant contribution to the thermodynamic 
properties to be calculated. Indeed, in usual approaches to the dilute Bose gas these terms 
are always neglected for this very reason. One aspect of the renormalization group cal- 
culation is that it is possible, and even rather straightforward, to include the three-body 



interaction in the calculations. As its influence will turn out to be extremely small, except 
in the critical region where it becomes somewhat larger, we will in the rest of this section 
omit this term for reasons of brevity. However, in Sec. Ill we perform some calculations 
including the three-body term to show its effect quantatively, and at that point introduce it 
again. 

Expanding the flelds in Fourier modes through 

and the complex conjugate expression for ?/'*(x, r), we can write the action in momentum 
space as 

1 1 

S[a*, a] = Y^-iTlUn + Ck - /iX ,,ak,„ + T^TI^ Y1 ^q'^k+q,n+m«k'-q.n'-m«k',n'ak,n • (4) 

k,n ^ripv ^ ^ 

n,n' ,m 

In this equation et = ^^k^/2m is the kinetic energy, Vq = J dx. V{'K)e~^'^'^ is the Fourier 
transform of the interaction potential, V is volume of the system and the bosonic Matsubara 
frequencies Un = 2Tm/h(3 reflect the periodicity of the flelds in imaginary time. 

The renormalization group equations now follow from repeatedly applying the renormal- 
ization group transformation to this action. It consists of three different stages [^,0 which 
we summarize here. The flrst step is to perform the functional integral in Eq. (1) for the 
most rapidly oscillating flelds. In other words, we eliminate the highest Fourier modes from 
the action in Eq. @). In practice, we therefore have to split the flelds in 'slow modes' and 
'fast modes', i.e. we deflne 

(x, r) = j^^^ Ek,n ak.ne*^''-"-""^) for < | k | < A - rfA 
V^>(x,r) = ^^^^Ek,„ak,ne*^''™"^ for A-rfA<|k|<A 
where A is the ultraviolet cutoff in the theory, and dA is the thickness of the shell in 
momentum space which is going to be integrated out. We will dwell on the issue of the 
ultraviolet cutoff shortly. With the deflnitions in Eq. (^, we can write the action as 

s[rM = So[r<,i^<] + So[r>,^>] + Si[r<,i^<,r>,i^>], (6) 



where So is quadratic in the fields, and Sj contains the remaining terms, and in particular 
the terms that couple the slow and fast modes. The partition function can now be written 

as 

Zgr = I ci[^*]c?[^<]exp|-i^o[^<,^<]} X J rfiV'* ]rf[V'>] exp { - ^ (5o[V^^ , ^>] + , , ^>] | 

= I d[r<:W<\^^p{-lso[r<,^<]} X cxp[-^As[r<,^<]] . (7) 

After performing the last integral, we obtain an expression for the partition function con- 
taining the new effective action S'[ijj^,ip^ = Sof^/'o V'<] + for the slow modes. 
The effect of this mode elimination is to change the value of the existing coupling constants 
and to introduce new vertices in the original action. The discussion of the rules on how we 
can obtain the equations describing this change of the coupling constants in general will be 
deferred to Sees. Ill and IV where we will also apply these rules to find the fiow equations 
for the Bose system. The two remaining stages of the renormalization group transformation 
consist of a rescaling of the momenta such that the new cutoff, which is A — dA, is restored 
to its initial value A, and a rescaling of the frequency and fields in the action such that 
there is no effect from this momentum rescaling on some appropriate terms in the quadratic 
part of the action. If we neglect the renormalizations from the first step, this procedure 
yields the so-called trivial scaling of the coupling constants and reveals the relevance of the 
various vertices as it tells whether the value of the corresponding coupling constant grows 
(termed relevant), remains the same (termed marginal), or shrinks (termed irrelevant) as 
we scale the ultraviolet cutoff. Having completed these final steps of the renormalization 
group transformation, we can read off the new coupling constants from the action. Thus, 
the renormalization group equations have two contributions. One is from the rescaling, the 
other is from the actual integrating out of high momenta from the action. 

Concerning the trivial scaling, there are two cases to distinguish. For the normal renor- 
malization group equations this scaling is found using the full quadratic action. The other 
situation occurs exactly at the critical temperature. At that point the correlation length ^ of 
the system is infinite, and the system looks the same on every length scale. Hence, the cou- 
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pling constants remain the same when performing the renormahzation group transformation 
and we are at a fixed point. To find this fixed point from the set of renormahzation group 
equations, we generaUy have to use a different scahng from the one found above, and this is 
the scahng at the critical temperature. It can easily be obtained by realizing that the corre- 
lation time Tc also diverges as we approach the critical temperature, behaving like Tc oc 
At temperatures such that Tc, we can neglect the time dependence in the action and 

the d/dr—term in the quadratic part of the action is unimportant. As a result, we then only 
need to consider the kinetic energy term to find the trivial scalings. This yields the scahng 
at the critical temperature. In renormahzation group studies, one normally restricts oneself 
to couplings which are relevant or marginal at the critical temperature, and this we will in 
first instance also do. However, it turns out that even the marginal coupling constant Uq 
from the three-body interaction term is quite irrelevant to the nonuniversal properties to be 
calculated. This is of course due to the fact that three-body interactions are unimportant 
for a dilute system. However, if one is interested in the universal properties of the phase 
transition, such as e.g. the critical exponent u connected to the divergence of the correlation 
length when approaching the critical temperature, marginal and also irrelevant couplings 
can have a considerable effect. We will explicitly encounter this fact in the next section. 
Also the position of the fixed point is changed when including irrelevant couplings. 

The renormahzation group equations describe the change of the coupling constants as 
we integrate out momentum shells from the action. Having arrived at the long wavelength 
effective action, we find whether or not we are in the condensed phase, and we can calculate 
universal properties connected to this phase transition. However, our aim is to use the renor- 
mahzation group method to obtain information on other properties as well. In particular, we 
want to calculate the equation of state, the superfiuid density, and the pressure of the gas. 
This can also straightforwardly be done by noting that e.g. for the total density we have 
n — Z^k,n(^k,n^k,n)- This equation can be cast into a differential equation describing the 
building up of the density as we proceed with the elimination of the fast Fourier modes. For 
that purpose, we have to use the right value of the chemical potential /i in each subsequent 
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momentum shell as found from the renormalization group equations. In the same fashion 
one can also determine e.g. the superfluid density and the thermodynamic potential. 

Before we can start with the derivation of the flow equations, we first have to pay some 
attention to the high momentum limit of the action. In principle, there is no real sharp 
ultraviolet cutoff A in our problem. However, the typical behavior of the Fourier transform 
of the two-body interaction potential, depicted in Fig. 1, is such that there is an effective 
ultraviolet cutoff around the momentum scale set by the scattering length a of this potential. 
Below this value, the Fourier transform is practically momentum independent and equal to 
Vq. As in the Bose systems considered here and realized experimentally the temperatures 
are so low that h/a ^ h/Afh, the particles in the gas reside in a momentum range well below 
this ultraviolet cutoff. Thus, we can represent the interaction potential by the momentum 
independent value Vq for momenta below a cutoff hA of 0{h/a), and zero for larger momenta. 

Modelling the potential as such implies that the nonuniversal properties we find from 
a renormalization group calculation will be sensitive to the specific value of the cutoff A 
taken in the calculations. However, at this point our knowledge about the microscopic de- 
tails of the Bose gas comes in to resolve this potential problem. In particular, we know 
that the two-body interaction potential has to renormalize to the two-body T— matrix 
T2^((k - k')/2 + q, (k - k')/2; h^(k - k')V"^) when we include all possible two-body scat- 
tering processes in the vacuum [^. The two-body T— matrix roughly has the same momen- 



tum dependence as V^, cf. Fig. 1, and is in particular constant and equal to Anah^ /m in the 
range of thermal momenta and energies. Thus, given an ultraviolet cutoff A we can fix the 
renormalization group equations by demanding that for the two-body problem, Vq indeed 
correctly renormalizes to Anah'^/m. Since this value is, due to the inequality a/Kth 1, 
already attained before we enter the thermal regime as we integrate out more and more mo- 
mentum shells from the action, this indeed leads to a correct description of the properties of 
the Bose gas which is independent of the ultraviolet cutoff A. Having eliminated the cutoff 
dependence, we are then in a position to determine also the nonuniversal properties of the 
dilute Bose gas. Furthermore, we can perform the calculation for any (positive) value of the 
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scattering length, thus being able to describe any atomic species with effectively repulsive 
s— wave scattering. The results we find are therefore relevant to the current experiments 
using ^^Rb and ^'^Na, but also to future experiments using atomic hydrogen or other atoms 
with a positive scattering length. 

We now turn to the application of the renormalization group method to the dilute Bose 
gas. First, we derive the renormalization group equations when the chemical potential is 
negative, starting from the action in Eq. (^. These equations however, do not describe the 
Bose condensed phase for then it is required that the chemical potential be positive. As 
a result, we have to rederive the renormalization group equations for that case and this is 
carried out in Sec. IV. The derivation is now much more involved as the space and time 
independent part of the effective action, i.e. — /i | ao,o P +Vq \ ao,o |^ /2hpV, has for /i > a 
Mexican hat shape and we first have to expand the action around the correct extremum by 
performing the shift ao,o aofl + V^ofi'PV and introducing the condensate density no- Only 
after that can we proceed to find the contributions to the renormalization of the vertices. 

III. THE SYMMETRY UNBROKEN PHASE 

In this section we first concentrate on the renormalization group equations valid for 
negative chemical potential. The reasons for this are threefold. The equations correctly 
describe the Bose gas in the uncondensed phase, and they offer an easy way to determine 
the influence of three-body effects on the quantities of interest. Moreover, it is best to start 
with this relatively simple set of equations because it illuminates most clearly our procedure 
to eliminate the cutoff dependence of the theory. 

A. The flow equations 

To calculate the change of the couplings after each step of the renormalization group 
transformation we can technically proceed in two different, but equivalent ways. The first 
and probably most familiar method is to expand the integrand in the partition function in 
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powers of S'/['?/'^, V'<, V'>5 V'>]5 ^ind then perform the integrals over the fast modes by evalu- 
ating the appropriate Feynman diagrams contributing to the renormalization or flow of the 
vertices of interest. The renormalization of the gradient and time-derivative terms of interest 
are found by performing a Taylor expansion in external momenta and frequencies of these 
diagrams. However, the associated couplings become more and more irrelevant for higher 
order terms in this expansion. The integration over the internal momenta of the diagrams 
is restricted to the afore-mentioned momentum shell, which can in principle be chosen to 
have any thickness dA. Taking it infinitesimally small, possible in the thermodynamic limit 
only, leads to a continuous renormalization group transformation described by a coupled set 
of first order differential equations. A much more important reason to take dA infinitesimal 
is that in this case the type of diagrams which contribute to the renormalization are the 
so-called one-loop diagrams only. This is a point which is not always explicit in the litera- 
ture, but it follows from the fact that each extra loop, i.e. each extra momentum integral, 



introduces an extra factor of dA ||30|, and therefore will vanish when dA is infinitesimal. 
Clearly, this is not so when the thickness dA has a finite value. Note furthermore that the 
above holds irrespective of the magnitude of the interactions, in contrast to what is some- 



times mentioned in the literature P^J2^. In the case of e.g. a large two-body interaction 



vertex, the infinitesimal thickness of the momentum shell still ascertaines the validity of 
a perturbation expansion with the use of one-loop diagrams only. This most notably im- 
plies that in this formulation the restrictions of the renormalization group method lie in the 
number of vertices considered, and not in the type of diagrams taken into account. The 
second way to obtain the renormalization group equations does not explicitly make use of 
Feynman diagrams, and is therefore very useful and efficient when the number of diagrams 
is large and/or the associated combinatorics is complicated. We here use the diagrammatic 
method as in the unbroken phase the number of diagrams is limited and the combinatorics 
is simple. We use the other method in Sec. IV. Moreover, using Feynman graphs also gives 
a transparent way to find the trivial scaling of the vertices in a somewhat different fashion 
than from the rescaling procedure described in Sec. II, which brings out more clearly the 
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physics of this procedure. 

As mentioned before, we include in first instance also the three-body interaction term 
written down in the previous section into our considerations, and take it, as the two-body 
interaction, to be momentum independent below the cutoff A. The one-loop diagrams to 
be calculated are depicted in Fig. 2. We find that in the thermodynamic hmit the total 
contribution to the chemical potential from integrating out an infinitesimal momentum shell 
in the Hartree and Fock diagrams is 

JA{i)-dA {2Try 

where A{1) — Ae~' denotes the radius of the shell in momentum space. The ladder and 
bubble diagrams renormalizing the two-body interaction potential give, together with the 
diagram containing the three-body term, 

,V„ ^ -Vi r + - AVi r ^,me. - ,m^, - „ + 1] + 

JA{i)-dA {27ry 2(ek-//) JA{i)-dA {27ry 

3^0 / 7^^(ek - 1^) , (9) 
JA{i)-dA {27ry 

and for the vertex Uq we find 

dUo = 8V^ / (3^N{e^ - /.) [N{e^ - /^) + 1] [2iV(ek - /^) + 1] + 

JA{l)-dA [ZTT)'^ 
rA{l) I 

/ JTT^J ^[1 + 2iV(ek + 2/3(ek - i^)N{eu - l^)[N{eu - /^) + 1]] - 

JA{l)-dA [ZTT)'^ (ek — /J') 

rm A 3(l + 2iVfa-,)) ^ 

J A[l)-dA {2'kY Ck-A* 
In these expressions iV(ek — jj) = l/(e^^^''~'^^ — 1) is the Bose-Einstein distribution function 
which result from the summation over the Matsubara frequencies 

To derive the renormalization group equations and the trivial scalings of the coupling 
constants we focus on the first of these equations renormalizing the chemical potential, and 
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cast it into a differential equation. Using | k |= A(/) = Ae~' and performing the angular 
integrals, we can write Eq. (||) as 

dfi = -2Vo-^ / N{eAe-^' - i^)e-^'dl . (11) 

/TT^ Jl 

Next, we remove all explicit /—dependencies from the Bose-Einstein distribution function 
by letting both the temperature and the chemical potential scale with exponent 2, i.e. we 
put T{1) = Te^' and fi{l) = fie'^K Hence, both temperature and chemical potential scale 
trivially with exponent 2 and we find for the differential equation describing the change of 
the chemical potential when integrating out a momentum shell 



^ = 2/i - —Voe-'N{eA - /i) . (12) 

Finally, we now also absorb the factor into Vq in order to remove the remaining explicit 
/—dependence. As a result, Vq scales trivially with exponent —1. By considering in Eq. (^ 
the term containing Uq, we can analogously show that this vertex trivially scales with an 
exponent —4. Note that to find the real physical quantities we should always remove the 
trivial scalings again. We see that the trivial scaling can, in a very simple way, be found 
from the one-loop expressions. Moreover, it shows that introducing the trivial scalings does 
not have an essential effect on the renormalization of the vertices; it is merely a rewriting of 
the differential equations. In the case of a negative chemical potential we eventually obtain 
the following coupled set of renormalization group equations for the coupling constants /i, Vq 
and Uo, 



da A^ 

dl TT 



2/^ - -Z^^oNieA - /i) , (13a) 



dVo_ A3 
~dr-~^'~ 2^^° 
3A3 
2^ 



^^UoN{eA-fi) , (13b) 



dU 

= -4t/o + [sP'N{eA - /i)[iV(eA - /i) + l][2N{eA - /i) + 1] + 
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[1 + 2iV(eA - /i) + 2/3(eA - fi)N{eA - /i)[iV(eA - /i) + 1]] 




3(1 + 2iV(eA - /i)) 



+ 18/?iV(eA-/i)[iV(eA-/i) + l] 



(13c) 



Ca - /Lt 



To argue that these equations are the only ones we need to consider, we still have 
to determine the trivial scaling at the critical temperature. This different scaling comes 



A^(eA ~ jj) = l//9(/)(eA — At(0) ^"^^ effectively at very high temperatures due to 

/?(/) = /?e~^'. To remove again all explicit /—dependencies after substituting this behav- 
ior, we clearly need a different trivial scaling of the vertices, and this is precisely the trivial 
scahng at the critical temperature, since putting A^(eA — /i) = l//3(eA — yu) is equivalent to 
neglecting the time-derivative in the action. In this manner we can straightforwardly show 
that the scaling of the chemical potential remains the same, i.e. we have /x(/) = /ie^', that the 
scaling of the two-body interaction becomes Voil) = Vbe' instead of Vb(0 = Vbe~', and that 
the three-body interaction does not scale, i.e. Uo{l) = Uq instead of Uo{l) = Uoe'^K From 
this we conclude that yU and Vq are relevant and Uq is marginal at the critical temperature. 
Furthermore, the four-body interaction is indeed irrelevant and therefore not included in the 
calculations. The coefficients of the gradient and time derivative terms in the quadratic part 
of the action are, like ?7o, marginal and would in principle also have to be included in the 
renormalization group equations. However, as we are in the regime a/Afh ^ 1, the interac- 
tions are expected to be independent of momentum and energy (but see below) and there is 
no renormalization of the d/dr and | V p terms. Moreover, we know from the e— expansion 
that the anomalous dimension rj, indicating the importance of the | V p renormalization at 
the critical temperature, is very small, namely rj = 0.02 [^. Therefore, the d/dr and | V p 
renormalizations will be neglected and Eq. (13) describes the renormalization of the vertices 
we will consider. 

To calculate the partition function of the gas, we need, next to the flow equations for 
the above quantities, also the correct boundary conditions. The first one for the chemical 
potential fi is just the bare value in the action Eq. (^. For the two-body interaction potential 



about because in the limit / 



oo the Bose-Einstein distribution function behaves as 
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Vo we need to be more careful, as this vertex has to correctly fix the renormalization group 
equations as described in the previous section. From Eq. ( |13b| ) we recognize that in a 
vacuum, i.e. A^(eA — /i) = 0, the renormalization of the interaction between two particles is 
governed by 

i[r-~^'~2^^'WZ^)- ^^^^ 

This is just the differential form of the Lippmann-Schwinger equation |^ for the two-body 



T— matrix at energy 2/i, i.e. T^'^(0, 0; 2/i). As the two-body T— matrix is energy independent 
for low energies, we can neglect /x and use T^^(0, 0; 2yu) ^ T^^(0, 0; 0) = ATrafi^ /m. As the 
solution of Eq. ([T^ ) is also practically independent of the chemical potential, we can there 
also neglect it. As a result, the requirement is now that, given an ultraviolet cutoff A, Vq 
flows for / — > oo to the value Anah^ /m. This can be ascertained by choosing the right initial 
condition for Vo, and more precisely we flnd from analytically solving Eq. ([1^) that 

Anah^ 1 , . 

Vo = = — (15) 

°^ ' m l-2aA/7r ^ ^ 

leads to the correct result. Note that we can describe different atomic species by only chang- 
ing the value of the scattering length a used in this equation. Finally, we in principle also 
need a boundary condition for the three-body interaction Uo- For this interaction we can, 
analogous to the scattering length a for the two-body interaction, introduce a length scale 
b, and again flx Uo{l = 0) such that the renormalization group calculation gives the correct 
result Uo{l = oo) = Anh'^b'^/m for elastic three-body scattering in a vacuum. However, in 
general not much is known about the microscopic details of the three-body interaction in a 
dilute Bose gas, and in particular about the value of the 'three-body scattering length' b. 
But, as can be expected from the fact that the three-body interactions are in the renormal- 
ization group sense irrelevant at large momenta, the results are practically insensitive to the 
boundary value of Uo used, and Uo{l = 0) is hardly of any importance. This is shown ex- 
plicitly in Sec. IIIB where we analyze the results from the renormalization group approach 
for ^^Na. Note, that taking Uo{l = 0) = is equivalent to assuming that three-particle 
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scattering is solely due to the sum of pair interactions. This is a standard approximation in 
atomic three-body calculations. 

Finally, to describe the dilute Bose gas we still need to derive expressions for the total 
and superfluid density, and the equation giving the thermodynamic potential Q, and thus 
the pressure p = —Q/V = —ui. To do so, we make use of the following well-known one-loop 
expressions for the density 

the superfluid density 

n^=n-nn = n- — — -/?ekA^(ek - /i)[A^(ek - /i) + 1] , (17) 

where n„ is the normal density given by the momentum-momentum correlation function, 
and the thermodynamic potential 

Casting these equations into a differential form by performing the integration shell by shell 
leads to 



dl 27r2 
dn^ dn h? 



dl dl Svr^ 

duj 1 A3 



/5eAiV(eA - yu)[iV(eA - ^^) + Ije'^' (19b) 
2/r^(l-e-^(^^-^))e-^^ (19c) 



dl P2n 

where the inverse temperature again scales as /?(/) = /3e~^' and the chemical potential is 



found from Eq. ( |13a| ) at each step of the integration. These equations describe the building 
up of these quantities as we integrate out more and more momentum shells from the ac- 
tion. For convenience, the explicit /—dependence is not removed from these equations, and 
they thus immediately yield the physical quantities. Note furthermore that these equations 
have no influence on the renormalization of /i, Vb, and Uq as they are decoupled from the 
renormalization group equations (13). 
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B. Analysis of the flow equations 



We start our analysis by first focussing on tlie critical properties of Eq. (13) and in par- 
ticular on the critical exponent v pertaining to the divergence of the correlation length, i.e. 
the correlation length behaves as ^ = I ^ Tcil^c when approaching the critical 
temperature. For that purpose we have to find the fixed point of the renormalization group 
equations, linearize the flow equations around this fixed point and identify the largest eigen- 
value A+ which is related to this critical exponent via v = 1/A+ ||^. We in first instance 
omit the three-body interaction, but in a subsequent calculation include it again to deter- 
mine the influence of this marginal vertex on the value of u. Only for the set {/i, Vo}, when 
Uo{l) = 0, do we perform the calculation of the fixed point explicitly. With the remarks 
made in Sec. IIIA we have that the fixed point is found from 

•^ = 2^-^V,J^ = (20a) 

yielding (/i*,'Ko*) = (eA/6, STr^e^/lSfcsTA^). From linearization of the Eqs. ( |20a|) and ( P0b|) 



around this value we find for the largest eigenvalue A+ = 1.878, implying that u = 0.532. 
Repeating the calculation including the equation for Uo{l), the fixed point is shifted and 
the critical exponent is found to be z/ = 0.613. Thus, we see that the marginal operator Uq 
has a rather large effect. Moreover, we can conclude from this result that also irrelevant 
coupling constants must have a considerable effect as it is known, from the e— expansion 
of the 0(2)— model and from measurements in ^He experiments, that the true critical 



exponent of the Bose gas should have the value u = 0.67. This discrepancy should be 
alleviated by including more and more irrelevant vertices. 

However, as we are in particular interested in the nonuniversal properties of the dilute 
Bose gas, we now turn to the influence of the three-body interaction term Uq on these 
quantities. The influence of this term will of course be largest when we start close to 



the critical chemical potential, because at the critical point in principle all fluctuations 
are of importance. Starting with a chemical potential near the critical value leads to a 
trajectory that almost flows into the flxed point of the renormalization group equations, 
and the momentum interval in which irrelevant vertices can have a signiflcant contribution 
to the flow, and also to the building up of the density, is then largest. The bare chemical 
potential yielding a flow into the flxed point is positive, and the physical chemical potential 
(i.e. with the trivial scaling removed) renormalizes to zero. For a bare chemical potential 
larger than this critical value the flow is no longer deflned as at some value of I we have 
that (eA — A*(0) becomes zero and the Bose-Einstein distribution function diverges. We will 
come back to this point later on and restrict ourselves here to the accessible regime, which 
physically implies that n < ric- 

The flrst aspect connected to the three-body interaction concerns the initial value prob- 
lem for Uq. Indeed, as alluded to before, an exphcit calculation shows that changing the 
boundary condition for Uq from to one corresponding with a 'three-body scattering length' 
b = 10a = 520ao, which is extremely large in general, changes the total density and the pres- 
sure in the system with less than 0.1 %. Thus, the results we obtain are practically insensitive 
to this boundary condition and henceforth we simply use b — a — 52ao. The second aspect 
we want to consider is the influence of the three-body interaction term itself on the outcome 
of the renormalization group flow. This we do by alternatively including and excluding this 
vertex. That is, we solve the set {//, Vq, C/q} and the set {//, VJ)} and compare the results 
we flnd. For that purpose, we plot in Fig. 3 the p — n~^— diagram near the critical density 
ric, where the influence of [/q is largest. From this flgure we see that the change in density 
and pressure is about 1% at maximum. Far away from the critical conditions, i.e. at large 
negative chemical potential, the system becomes more and more dilute, and the influence 
of the Uq term vanishes, consistent with expectations. In principle, we could choose to 
maintain the three-body interaction term in the renormalization group equations. However, 
as its effect is very small we will from now on neglect Uq altogether. Thus, the dilute Bose 
gas, and more in particular the nonuniversal properties we are interested in, is accurately 
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described by only following the renormalization of the chemical potential and the two-body 
interaction. Having concluded this, we restrict ourselves from now on to the coupled Eqs. 
( |13a| ) and ( [I3E| ). 

However, before we analyze some physical implications of these equations, we want to 
remark that the dependence on the ultraviolet cutoff A is indeed eliminated from the theory. 
The influence on e.g. the density can be shown to be completely absent, of course with 
the limitations that should be larger than the thermal momentum h/Kth and that Vq 
is properly renormalized to Aitah^ /m when we enter the thermal regime. We are going to 
compare the results from the renormalization group calculation with known results for the 
weakly-interacting Bose gas as found from the many-body T— matrix theory [p!5| , p2| . Far 



from the critical temperature we expect the results of the mean-field and renormalization 
group calculations to be identical. However, close to the critical temperature the renor- 
malization group method will clearly deviate from mean-field theory. We will here focus 
on the behavior of the effective two-body interaction, and defer the discussions concerning 
the equation of state and other thermodynamic quantities to Sec. IV when we are able to 
describe also the condensed phase of the gas. 

In the many-body T— matrix theory, the chemical potential is renormalized to ^' = 
jj, — 2nT^^(0, 0;0) ^2[|. This is one of the results of including all two-body scattering 
processes, but also incorporating the effect of the medium on the scattering. Including the 



latter effect on the collisions in the gas, which was first carried out explicitly in Refs. |q], |15 



and [0 , is an important step forward in the correct mean-field treatment of the dilute Bose 
gas, since including many-body effects causes the effective interaction to go to zero at the 
critical temperature. This resolves a number of problems found in previous approaches using 



just the two-body scattering length a [|T6|-|T^. With our renormalization group approach 
we can corroborate this result and even go somewhat further than that. In our previous 
work we included only the many-body effects coming from the ladder diagrams. However, 
a class of diagrams that in principle also affects the two-body interaction are the bubble 
diagrams. With our current set of renormalization group equations we can precisely study 
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the effect of these bubble diagrams on the effective interaction. This is straightforward 
because we can pinpoint the ladder and bubble contributions in the equation describing 
the renormalization of the interaction. By alternatively including and excluding the bubble 
diagrams and then solving the renormalization group equations we can study the relative 
importance of the bubble diagrams on the effective 'many-body' scattering length a^^^ . To 
avoid any confusion we will adopt the following notation for the scattering length. The 
normal 'bare' two-body scattering length as found e.g. from analysis of the appropriate 
association spectra is denoted by a, as usual. The effective scattering length a^^^ includes 
also effects of the medium on two-body scattering, and therefore depends on the specific 
approximation used to calculate this effect. Here we consider two such approximations and 
denote the corresponding scattering lengths by a^^ and a^'^ . In the first case it is the result 
of a many-body T— matrix calculation and is defined through T*^^(0, 0, 0; 0) = Aira^^Ji^/m. 
In the second case it is the result of a renormalization group calculation and is defined 
through Vo{l — oo) = Ana^h'^ /m, irrespective of the fact if bubbles are or are not included 
in this calculation. In Fig. 4 we depict the ratio of the effective scattering length resulting 
from the renormalization group approach to the simple two-body scattering length a when 
we include and exclude the bubble diagrams, as a function of T/T^ and at a density of 1.5 10^^ 
cm~^. We conclude that the effect of the bubbles can be rather large and is in particular 
important near the critical temperature. Further away from the critical temperature the 
importance of the bubbles rapidly decreases. This imphes that calculating the influence of 
many-body effects by means of the ladder diagrams (i.e. doing the full many-body T— matrix 
calculation) can quantatively give a poor estimate for a^^^ . However, qualitatively there 
is clearly good agreement as we find from the renormalization group calculation that the 
effective scattering length indeed goes to zero at the critical temperature. 

As already mentioned before, a bare chemical potential larger than the critical one, 
which is positive and yields a flow into the fixed point, corresponds to an inaccesible density 
regime. Therefore, we are not able to penetrate the region with a Bose condensate and 
cannot describe the Bose gas below the critical temperature with the renormalization group 
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equations in Eq. (13). Moreover, the critical exponent v found from this set is not a very good 
approximation to the true value v = 0.67, even if we would include three-body interactions. 
These aspects are intimately related and due to the fact that in the case of a positive chemical 
potential, the k = part of the effective action has a Mexican hat shape. Therefore, we 
must explicitly break the symmetry and introduce the condensate density into the action 
by expanding the action around its minimum, and not around (ip) = as was done in this 
section. For a negative chemical potential, the above approach is of course correct as then 
{ip) is equal to zero. 



IV. THE SYMMETRY BROKEN PHASE 

Breaking the symmetry allows us to describe the dilute Bose gas below the critical tem- 
perature. Moreover, also above the critical temperature we find considerable improvement. 
This is due to the fact that a positive bare chemical potential can renormalize to negative 
values. Thus, starting out in the broken phase, the fluctuations can restore the symmetry 
and we end up in the unbroken phase above the critical temperature. The new set of renor- 
malization group equations explicitly takes this broken symmetry into account and therefore 
gives a much better description than the one resulting from the set of equations in Eq. (13). 



A. The flow equations 

In order to expand the action around the correct extremum we have to perform the shift 



flo,o oo,o + y/noh(3V which introduces the condensate density uq. Substituting this in the 
action Eq. (^ leads to the familiar expression |^| 

S[a, a*] = -h(3V{fino - ^nlVo) + (-/ii/^ + nQy/n^VQ){al Q + ao,o) 
+ + Ck - + 2noV^o)ak,„ak,n + -UoVq ^( 

'-^k,n'-^— k,— n ~^ '^k,n'^— k,— n) 

k,n k,n 



+ y fl(3V ^ ^o(%,m'^k-q,n-m'^k,n + '^k+q,n+m^q,m'^k,n) 



n,m 
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(21) 

k,k',q 

n,n',m 

for the action. The magnitude of the condensate is determined by ehminating the hnear 
term from the action. In first instance we thus find no — /^/Vq. As a result, we can write 
the action as 

S[a, a*] = -h(3VuQ + ^(-ifiu;^ + et + rii)ak „ak,„ + ^Fia X]("k,n«-k -n + Ok.nO-k -n) 

k,n k,n 



/ fc/3T/ X/(^q,m^k-q,n-m^k,n + ^k+q,n+m^q,mCtk,n) 



vWk„ 



+ 2 fll3V ^ ^k+q,n+m^k'-q,n'-m^k',n'a.k,n , (22) 
n,n',m 

introducing a;o = nQVo/2 as the fowest order approximation to the thermodynamic potential 
density, and defining the vertices Fn — ji — ^Eii(0; 0), = ?iEi2(0;0) and r3 for which 
we have in lowest order that Fn = = uoVq — n and Fa = y/noVo — \^ fjVo- These are 
the boundary conditions for the flow equations we derive next. 

Deriving the renormalization group equations using Feynman diagrams is in this case 
more involved than in the unbroken phase as the number of vertices is larger, but more so 
because we now also can have anomalous propagators {dk^n'^-k-n) (ok.nO-k.-n) in these 
diagrams. The number of diagrams is therefore much larger and the combinatorics is more 
complicated. Therefore, we will here use a different method to obtain the renormalization 
group equations which does not explicitly make use of Feynman diagrams, and is therefore 
very useful and efficient in this case. It relies on the fact that in Eq. (7) 

I d[r>w>]ew[-l{So[r>,^>] + Si[r<,i^<,r>,^>])} = ew{^^^ , (23) 

where the trace is over Fourier modes in the shell dA only and G> is the Greens function 
for these fast modes. Taking the shell infinitesimal again, we can simply calculate ln(— G>^) 
because we then only need the part of the total action which is quadratic in the fast Fourier 
modes (yielding all one- loop contributions). The coefficients of this quadratic part also 
contain the slow modes because the interaction term Si[ip'^,ip<:,ip!^,ip>] couples the slow 
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and fast modes. By simply Taylor expanding ln(— G^^) we straightforwardly find the new 
effective action for the slow modes and the renormalization group equations of any vertex 
we would like to consider. 

Thus, we split the fields in slow modes and fast modes as in Eq. (||) and find that the part 
from the action which is quadratic in the fast modes and only leads to the renormalization 
of coupling constants reads 

^^'^[V^>,V^>] = E'(-«^^n + ek + rn + 2r3(V^* +V^<) + 2Vo I r)a^,„ak,„ 

k,n 



+ (^ri2 + r3^< + ^Vo^^)E' 



k,n^-k,-n 
k,n 

+ (^ri2 + Tsr< + ^Vo^*')E «k,na_k,-„ , (24) 

k,n 

where the prime denotes that the sum over momenta is restricted to an infinitesimal momen- 
tum shell dA at the cutoff. Evaluating the functional integral over these fast modes leads to 
adding Tr(ln(— G^^)) to the action for the slow fields, and thus changes the vertices. This 
quantity is most easily evaluated by performing a Bogoliubov transformation to diagonalize 
Eq. g, and we find that 



27r2 



^ ln(l - e-^^^) + ^(^A - (eA + En + 2E3(^* + ^<) + 2Vo | H) 



dK , 



where the second term originates from the diagonalization procedure, and the 'dispersion' 
E\ is found from 

El = (eA + Tn + 2T:,{r< + V^<) + 2V^o | ^< ?f - 

(Ei2 + E3^< + ^Vb^^)(Ei2 + T,r< + \vor<) ■ (26) 

In zeroth order in ip^ and tpt we retrieve the well-known Bogoliubov dispersion hujjy = 



(eA + Eii)2 — Efg, equal to hu/^ = + 2EiieA and thus gapless at Z = due to the 



equality Eii(/ = 0) = Ei2(/ = 0). This corresponds to the Hugenholtz-Pines theorem p5 
Performing a Taylor expansion in terms of the slow modes we find the new action. Thus, 
integrating out a momentum shell renormalizes the existing vertices, but also generates new 
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terms in the action, and in particular a linear term. To eliminate this term and remain 
in the minimum of the action we again have to perform a small shift in ao,o- This implies 
that the magnitude of the condensate changes as we are integrating out momentum shells 
and we also have a flow equation for the condensate density. Since we omit three-body 
interactions containing six fields, we also have to neglect terms containing five fields since 
they correspond, together with the condensate field, with a three-body term. As a result, the 
action remains of the form written down in Eq. (^) and the renormalization group equations 
can now be obtained for all the vertices of interest. However, due to the ?7(1)— symmetry 
of our problem, we can relate some of these vertices and thereby limit the number of flow 
equations we actually need to describe the complete renormalization of the action in Eq. 
(p^). As this f/(l)— symmetry cannot be broken during the process of renormalization, the 
action can, at any time, be recast in the explicitly U{1) symmetric form 



^<(x,r) + iVo(0 |^<(x,r) \' 



From this it is then first of all easy to see that the Hugenholtz-Pines theorem holds, 
implying that rii(Z) = ri2(/) = (See appendix A for an exphcit derivation of this 

important relation.) Next, also due to the neglect of three-body interactions, we have 



that r3(/) = yno(/)Vo(0 and Vo(/) = ri2(/)/Tio(/). Thus, we only need to know the fiow 
equations for e.g. nQ{l) and ri2(/) and then the other renormalization group equations can 
be determined. In the following we can therefore restrict ourselves to the renormalization 
of the linear term and the term proportional to {ip*^ + ipl^ as these determine the change 
of the condensate density and the anomalous selfenergy respectively. Note that the 
equality rii(/) = ri2(/) ensures that the dispersion Tiuii^ is gapless at any point during 
renormalization, as it should. 

To arrive at the fiow equations for no(/) and ri2(/), we first of all need the linear term 
dVQ^aQ Q + ao,o) that is generated by integrating out a momentum shell. We find that 
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Analogously we find a change dT^2 the anomalous selfenergy. (See appendix A for details.) 
However, before we know the full renormalization of this vertex we have to determine the 
shift needed to eliminate the linear term that is generated. Substituting a^ Q — >■ ao,o + 
y/WJW d{y/n^) and retaining only the term linear in d{y/n^) we find that 

J- 11 + J- 12 

which influences the renormalization of because we have for the total renormalization of 
the anomalous selfenergy 



dTu = dT';!J + 2T^d{^) . (30) 

due to this shift. Using the above mentioned relations between the vertices implicate from 
the f/(l)— symmetry we then make contact with the renormalization group equations for the 
unbroken phase by determining the flow equations for the chemical potential fi = and 
the two-body interaction Vq = Tu/no. After some algebra we ultimately find 



df^ o T/ 



^'''\lTV^\ 2N{h.,) + 1) - 1 + ^^^''Y^' mnu^.mnu,) + 1] 

ZlX ft CJ^ 



dVo A3 



dl 2ti 



2 

2 "^O 



^^^{2N{hut:) + 1) + ^^^±^pN{nu^)[N{hu^) + 1] 



(31b) 



and the condensate density follows from no(/) = /i(Z)/Vo(/). Breaking the symmetry is 
irrelevant to the trivial scaling of the vertices, and thus these are identical to what was 
found in the previous section. Comparing these flow equations to Eqs. (|13a|) and ( |13b|) 



omitting f/o, we see that both sets coincide when /i is taken equal to zero. Thus, the 
renormalization group equations for negative and positive chemical potential yield a flow 
which is continuously differentiable, also at yU = 0. 

We now know the renormalization of the chemical potential and the two-body interaction. 
Finally, we again have to find the equations describing the building up of the total and 
superfiuid densities and the thermodynamic potential as we are integrating out momentum 



(31a) 
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shells. Analogous to the procedure followed in the previous section we find for the total 
density, being the sum of the condensate density and the above condensate density, 

'^'^ (7^i^Nihuj^) + l)-l)e''^ (32) 



dl 27r2 \2niUA 2, 
with the boundary condition n(l = 0) = no(/ = 0) = [lil = 0)/Vo(/ = 0). The superfiuid 
density follows from 

^ = f - ^.pe^NihuA) [N{huA) + l]e-3' , (33) 
and the thermodynamic potential from 

^ = 2^ (^^"^' - e-/^^-)e-^' + lihu, -e,- ;.)e-3') , (34) 
(cf. Eq. (|25|)) with the boundary condition u{l = 0) = Uq = 

B. The critical temperature of BEC 

With the boundary conditions mentioned above we can again numerically integrate the 
renormalization group equations. For a fixed temperature, we vary the value of the (posi- 
tive) bare chemical potential, and calculate e.g. density and pressure. The physical effective 
chemical potential, i.e. with the trivial scaling removed, decreases as we perform the inte- 
gration, and depending on the starting value remains positive, renormalizes exactly to zero, 
or becomes negative at some value of the integration parameter /. In the first case we start 
out and remain in the broken phase and are below the critical temperature of Bose-Einstein 
condensation, i.e. we have a finite condensate density. In the second case, the condensate 
density rio(0 — A*(0/^o(0 renormalizes exactly to zero for / oo, and we are at the critical 
conditions for Bose-Einstein condensation. In the third case we started out in the broken 
phase, but the fiuctuations restore the symmetry. At the value of I for which the chemical 
potential becomes negative we have to continue the integration with the set Eqs. ( |13a| ) and 



and we are thus in the uncondensed phase. Hence, to be able to describe the dilute 



Bose gas above, but not too far from the critical temperature, we need the renormalization 
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group equations for both signs of fi. In Fig. 5 we depict the trajectories resulting from the 
integration of Eqs. ( |31a| ) and ( pib| ). 



It is evident from this figure that the critical properties of the Bose gas are determined 
by this set of equations. By linearizing the flow equations around the fixed point we can 
identify the largest eigenvalue A+, and determine the critical exponent u following from this 
set. We find u = 0.685, which gives a much better approximation to the critical exponent 
than the renormahzation group equations studied in Sec. Ill and is to be compared with 



the value u = 0.67 found from the e— expansion of the 0(2) model and measured in 
^He experiments. The agreement is surprisingly good, and together with the fact that we 
explicitly showed that three-body effects are negligible, this indicates that we are indeed 
accurately describing the Bose gas with the derived renormahzation group equations, also 
in the critical region. The cause of this good agreement is that although we only consider 
the renormahzation of fi and Vq, the type of scattering processes in terms of real (bare) 
particles we are actually taking into account are very elaborate. The propagator for the 
the Bogoliubov quasiparticles is namely the result of dressing the original bare propagator 
Ti/{iTiuJn — ca + /u) with interactions with the condensate as we use the terms 2noVoak n'^k,n 
and ^^o^o('2k n'^-k -n + 'j^k,nO-k,-n)/2 in the zeroth order quadratic part of the action. Thus, 
the diagrams we calculate actually contain an infinite number of scattering processes with 
the condensate. Therefore, we are describing the system much better than in Sec. Ill already 
at this level of renormahzation group. 

The first nonuniversal property we concentrate on is the change in the critical tempera- 
ture of Bose-Einstein condensation due to interaction effects. This result is presented also 



elsewhere [Q, but we will recapitulate it here because of its experimental significance. At 



fixed temperatures we vary the (bare) chemical potential to find the trajectories flowing into 
the fixed point. Using Eq. (p^), this yields the critical densities for these specific temper- 
atures and gives us the ric — T relation at which Bose-Einstein condensation occurs. We 
repeat this for different values of the scattering length to obtain the dependence of the 
critical temperature on the strength of the interaction. In Fig. 6 we show the degeneracy 
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parameter UcA^f^ found from the renormalization group calculation as a function of a/Ath- 
As seen from this figure we conclude that the critical temperature is raised with respect 
to the ideal gas value. This is in qualitative agreement with the recent experiments [jl],0] 
and also preliminary Quantum Monte Carlo calculations seem to confirm this result |^ 



An indication of an upward shift was also found some time ago by one of us studying the 
nucleation of Bose-Einstein condensation p|. From our calculations we predict that for the 
^^Rb and ^^Na experiments the critical temperature can be raised with as much as 10 %, 
which appears to be a very promising result because one might expect that an effect of this 
magnitude can very well be measured in future, more accurate, experiments. It is important 
to note that this shift in UcA^f^ can be observed in magnetically trapped atomic gases if one 
directly measures the density in the center of the trap at the critical temperature. One 
should in particular not measure the total number of particles, because this involves the 
density profile in the trap and due to the repulsive nature of the interactions thus tends to 
obscure the effect [ p^ . 

The reason for a higher critical temperature, or more precisely, a lower critical density 
is the following. The effective chemical potential renormalizes from a positive initial value 
exactly to zero. Consequently, we have the Bogoliubov dispersion in the equation for the 
density, and this depresses the occupation of the non-zero momentum states relative to the 
ideal gas case, where we would just have eA in the Bose-Einstein distribution function. The 
magnitude of the effect is related to the behavior of the chemical potential when renormal- 
izing to the fixed point value fi* = eA. Suppose we effectively have ^{l) = ae\ independent 
of /, with some positive a smaller than 1. We can then translate the differential equation for 
the density into an ordinary integral over momentum space following the inverse procedure 
from which we found the flow equations in Sec. III. Doing so, we find that the Bogoliubov 



dispersion effectively becomes hu]^ = vT+~2aek. Thus, this essentially boils down to a 
mass renormalization, and the change in the critical temperature is directly related to the 
magnitude of this renormalization. Starting with the equation describing the building up of 
the above condensate density (cf. Eq. (|A.6a|) from appendix A) one can easily show that we 
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approximately have UcA^f^ = (1 + a)g3/2{l)/{l + 2a). For an increase of 10 % in the critical 
temperature, we find from this result that a must be equal to 0.2. Indeed, this turns out to 
be the typical value of a in the thermal regime where the contributions to the density are 
largest. 



C. The region naAf^^ <^ 1 

To obtain also information on the properties of the Bose gas below the critical tem- 
perature we have to start with a chemical potential which always remains positive under 
renormalization. We are then always in the condensed phase and have to use the renor- 
malization group equations ( pla| ) and ( |31b|) . Written as ordinary one- loop integrals, these 



equations contain infrared divergencies and a straightforward perturbative analysis is not 
possible. This is a well known problem in the theory of the interacting Bose gas p8| , |39| . 
However, doing the calculation by means of the renormalization group approach, this prob- 
lem is in principle resolved due to the resummation which is automatically performed. The 
presence of infrared divergencies causes both the physical chemical potential and two-body 
interaction in the renormalization group approach to renormalize to zero, instead of becom- 
ing infinite as they would in a regular one-loop calculation. Indeed, one can show from Eqs. 
( plap and (|31b| ) that for / ^ oo we have oc e' and Vq oc e~^'. These scalings are different 
from the trivial scalings fi oc e^' and Vq oc e~', and are therefore termed anomalous. Phys- 
ically, the anomalous scaling implies an effective energy and momentum dependence of the 
coupling constants, as we show explicitly in appendix B. In particular, the anomalous scaling 
we find leads to the result that the chemical potential, and thus the normal and anomalous 
selfenergies, as well as the two-body interaction, behave linearly with k for low momenta, i.e. 
our renormalization group approach reveals that /iEii(k;0) oc k and /;,Si2(k;0) oc k. This 



behavior, implying that hTiii(0; 0) = /iSi2(0; 0) = 0, is an exact result for low momenta pS 
which we explicitly recover here. The consequences for the application of renormalization 
group are however twofold. First, it implies that the Bogoliubov dispersion hu^ does not 
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possess a sound mode anymore since the low momentum behavior is not hnear, but instead 



and ( |31b| ) are insufficient to describe the sound mode. Indeed, following the argument in 
Sec. Ill we neglected the renormalization of the (marginal) d/dr and | V ^ terms since the 
interactions were anticipated to be momentum and energy independent. However, below 
the critical temperature we see that the momentum and energy dependence of the selfener- 
gies does become important, and should therefore be included in the renormalization group 
calculation by Taylor expanding the selfenergies in terms of the external momentum and 
frequency. Put differently, the anomalous dimension rj is no longer small below the criti- 
cal temperature. The extra renormalization group equations obtained in this manner may 
change the particular anomalous scaling found above, but since ^Sii(0; 0) = ^Si2(0; 0) = 
is an exact result we still expect to have anomalous scaling of the coupling constants. The 
effect of the extra renormalization group equations will be to change the dispersion relation 
in such a way that the linear sound mode is recovered. We will not pursue this calculation 
here, but postpone it to future work. 

The reason for the anomalous scaling and consequently the disappearance of the sound 
mode, is caused by the infrared divergence in the one-loop expressions for the selfenergies 
and two-body interaction. This can be traced back to the behavior of the coherence factors 

and fk of the Bogoliubov transformation diagonalizing the quadratic part of the action. 
The Bogoliubov transformation is given by 



we have /itUk oc k^^"^. This is of course an incorrect result I^H] indicating that Eqs. ( |31a] ) 



-k,-n' 



(35a) 




Wk6k,n - 



^^k&-k, 



(35b) 



where 6k n 6k,n are the Fourier components of the fields creating respectively annihilating 



a Bogoliubov quasiparticle. The coherence factors are given by 




(36a) 
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(36b) 



and satisfy the requirement that ul- — v^. = 1 because the Bogohubov transformation is 
unitary. As the Bogohubov dispersion is given by ^tUk = + 2/iek, it is easy to see that 
Mk ^ 1 and f k — ^ for large momenta, and both factors behave as 1/ \/k as k tends to 
zero. The region where the crossover occurs between these two regimes is determined by the 
parameter naA^^, found from comparing the chemical potential with the thermal energy, 
and also determining the crossover from quadratic to linear behavior of the Bogohubov 
dispersion. When naK^^ ^ 1 the linear regime of the dispersion is extremely important in 
determining the properties of the Bose gas. Conversely, when nah^^ <^ 1, the Bogohubov 
dispersion can be approximated by the normal dispersion ek + /i, except for a very small 
region at low momenta, which is then not very important and can for most practical purposes 
be neglected. In the following, we therefore concentrate on the regime naAff^ <^ 1, covering 
a large temperature interval below the critical temperature, and use Wk = 1 and f k = 
in contributions to the renormalization of the vertices that contain infrared divergencies. 
Otherwise, we will use the full expression Eq. (|36| ) . This procedure turns out to be necessary 
to ensure the Hugenholtz-Pines theorem to be satisfied, as is shown in appendix A. 

In the method used in the previous section it is not possible to pinpoint the Bogohubov 
coherence factors at any stage of the derivation of the renormalization group equations. 
Therefore, we really have to go through the calculation of the Feynman diagrams relevant to 
the various vertices we are interested in. In our case these are the linear term of the action 
and the anomalous selfenergy. The one-loop diagrams of interest are depicted in Fig. 7 and 
contain also the anomalous propagators (Okn'^-k.-n) (flk.nCt-k.-n)- Using the Bogohubov 
transformation, these diagrams are straightforward to calculate, and we find that for the 
renormalization of the condensate density both diagrams for Fq contribute, but that for the 
renormalization of the anomalous selfenergy only diagrams A and B give contributions that 
do not contain infrared divergencies. The explicit calculation of these diagrams is presented 
in appendix A. Using the expressions obtained there, we find that the flow equations become 
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^ = 2/i - ^Vo {^-^j^{2N{hu^) + 1) - 1 + AP^^N{huJ^)[N{nuJ^) + 1]) , (37a) 

^ - {^^^ - M3mr^.mn^.^ . U) . (37b) 



where we again made use of the relations jj,{l) = Tii{l) = T^i^l) = ^no{l)T3{l) = "^0(0^(0 
due to [/(I)— symmetry. Note that also this set coincides with Eqs. ( |13a|) and ( |13b|) when 



the chemical potential is equal to zero, so the flow is also in this case everywhere continuous 
and continuously different iable. Moreover, the equations for the density, superfluid density 
and thermodynamic potential do not contain an infrared divergency, so the flow equations 
for these quantities remain the same and are given by Eqs. (|3^), (|33D and (p^D, respectively. 



D. Analysis of the flow equations 

At this point we can describe any point in the phase diagram of the dilute Bose gas. For 
negative chemical potential we have Eqs. ( |13a| ) and ( |13b| ), for positive chemical potential 



we must use Eqs. ( |37aD and ( |37b| ). Moreover, we have to combine both sets when we 
are not too far below the critical density when the chemical potential changes sign during 
application of the renormalization group transformation. Using Eq. (37), we find that the 
critical temperature of Bose-Einstein condensation changes with less than 0.1 % compared 
to the more accurate result found using Eq. (31). Therefore, this shows that using the set 
not containing the infrared divergencies essentially leads to the same results, implying that 
it is indeed correct to neglect the linear part of the Bogoliubov dispersion. Having proven 
this explicitly, we will in the following present the results of the renormalization group 
calculation of the effective two-body interaction, i.e. the many-body scattering length a*^-^-^, 
the condensate and superfluid densities and the p — diagram below as well as above the 
critical temperature, and compare them in all cases with the results from the many-body 
T— matrix calculation. 

We start with the scattering length a*^-^-^. Above the critical temperature it is straight- 
forward to take into account only the contributions of the ladder diagrams, or to include 
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the effect of the bubble diagrams as well. Below the critical temperature there are no clear 
ladder and bubble diagrams, so the various contributions are in principle not clearly associ- 
ated with ladders or bubbles. However, comparing the differential equations for positive and 
negative fi we can conclude that the first nontrivial term on the righthand side of Eq. ( |37b| ) is 
effectively related to a ladder contribution, and that the second nontrivial term is effectively 
related to a bubble diagram. We depict in Fig. 8 the scattering length a^-^-^, normalized to 
the two-body scattering length a = 52ao for ^^Na, when including and excluding the bubble 
contributions {a^'^ /a), and also the result found from the many-body T— matrix calculation 
{a^^ /a) [0. As our approach is restricted to the regime naAfi^ ^ 1, we present no results 
for nah^^ larger than one. The renormalization of the effective two-body interaction to zero 
at the critical temperature, a result already found in the many-body T— matrix approach, 
turns out to be correct. Indeed, this can be easily understood from the renormalization 
group equations. A fixed point (//*, Vq*) is present in the set Vb(/)e^'}, which means 

that for / oo the physical two-body interaction, i.e. with the trivial scaling removed, 
behaves as VQ*e~' and thus goes to zero. The depression in the scattering length around the 
critical temperature occurs in a fairly large temperature interval. Note furthermore that 
applying the renormahzation group equations ( |31a| ) and ( |31b| ) would lead to a^'^ = ev- 



erywhere below the critical temperature due to the anomalous scaling found from this set. 
This property is actually expected to hold true even in a more elaborate renormalization 
group calculation and is an issue worthwile studying as it may have important consequences 
for e.g. the exact form of the Gross-Pitaevskii equation describing the condensate. 

Next, we turn to the equation of state. For the ideal Bose gas we have in general that 

nAl = noAl + 93/2(0 (38) 
where gn{z) is a Bose-function defined as 

and ( = e^'^ the fugacity. In the many-body T— matrix theory above the critical temperature 
we have that the chemical potential is renormalized to fi' = fi—2nT'^^{0, 0; 0) |3^. Below the 
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critical temperature the dispersion changes to the Bogohubov dispersion and the equation 
for the density is in essence given by Eq. (0), recast in the form of an ordinary one- loop 
integral. In Fig. 9 we plot ^A^^^ for the ideal Bose gas, following from the many-body 
T— matrix calculation, and from the renormalization group calculation. In the latter case, 
we numerically integrate the flow equations for a fixed temperature, and vary the value of 
the bare chemical potential. This changes the total density in the system, and therefore 
the value of nA^^. For the ideal Bose gas, positive values of the chemical potential are not 
allowed. With the repulsive interactions taken into account, a positive value is possible 
and the fugacity can be larger than one. The renormalization group calculation yields the 
same density for a chemical potential slightly above and slightly below the critical chemical 
potential. This double valuedness in the density occurs in an extremely small region around 
the critical density and the effect is smaller than 1 promille in the case of ^^Na. For ^H, with 
a scattering length a = 1.34ao, the effect is even far below the promille level. However, it 
is important to note that this effect is much smaller than the change in density we obtain 
when we include three-body interactions, being approximately 1 % (cf. Sec. III). Therefore, 
the double valuedness we find cannot be trusted physically and should be neglected. It 
presumably disappears when we include three-body effects or extend the renormalization 
group calculation otherwise. 

When the bare chemical potential decreases the gas becomes more and more dilute and 
the influence of the interactions can better and better be accounted for using mean-field 
theory. This fact is actually evident from the renormalization group equations in Eq. (13). 
When the chemical potential is large and negative, the Bose-factors A^(eA — /i) are strongly 
depressed and there is hardly any many-body effect on the renormalization of the two- 
body interaction so we will find that it is just renormalized to T^^(0,0;0) = Anah^/m. 
Indeed, the same is true in the many-body T— matrix calculation |jl5[, since we have that 
T*^-^(0, 0, 0; 0) ~ T^^(0, 0;0) when the system is extremely dilute, i.e. there is no effect 
of the medium in this regime. In addition, the differential equations for /i and Vq are now 
almost decoupled, and consequently the chemical potential will renormalize approximately 
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to yU — 2nT^^(0, 0; 0), which is only a small effect in this limit. Therefore, the renormalization 
group equations can practically be recast in the regular one-loop expressions one encounters 
in the many-body T— matrix theory and the results we find in this regime are approximately 
the same. 

In Figs. 10(a) and 10(b) we depict the condensate and superfluid fractions as a function 
of temperature for a density of 1.5 10^^ cm~'^ both from the renormalization group and 
many-body T— matrix calculation. Again, not too close to the critical temperature we 
have good agreement. Note that the superfluid density has not yet become zero at the 
temperature where the condensate density vanishes. To explain this aspect we will focus 
on the equation for the superfluid density in the unbroken phase. When we interpret the 
fact that the chemical potential is renormalized as we integrate out momentum shells as 
a chemical potential depending on momentum, we can express the superfluid density as a 
regular integral over momentum space, i.e. 



Writing fj,{k) = /xq + Sfi{k) and performing a Taylor expansion, we find to first order in 6fi 



where the lowest order terms drop out. This is due to the fact that these terms yield the 
expression for the superfluid density in an ideal Bose gas above the critical temperature 
which can be shown to be exactly equal to zero. It is clear from Eq. (|4l|) that the superfluid 
density will in general not be equal to zero. The behavior of 6^{k) would have to be very 
special to give a superfluid density exactly equal to zero. It is therefore not surprising that 
we find from our limited set of renormalization group equations a different temperature 
for which the condensate and superfluid densities vanish. For the situation depicted in 
Fig. 10 we have however that AT/Tc is only about 8 10^'^. Extending the renormalization 
group calculation would in principle lead to a superfluid density which vanishes at the same 
temperature as where the condensate density becomes zero. 




(40) 



that 



r (i k 2 

J (2^^^^"^ " /^o)[iV(ek - /io) + 1][1 - 3/9ek(2iV(ek - /xo) + 1)] Sfi{k), (41) 
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Finally, we also present the pressure of the dilute Bose g function of the inverse 

density. Using Eq. (0), we depict this behavior in Fig. 11 for a ^^Na gas at 1 fiK, together 
with the result of the many-body T— matrix calculation. Although the critical densities are 
different, we find a fairly good agreement between the two curves and the difference between 
the renormalization group calculation and the many-body T— matrix calculation is small 
when it concerns the pressure of the gas. However, concerning the other nonuniversal prop- 
erties discussed in this article, the difference can be substantial near the critical temperature 
as we have seen. 



V. CONCLUSIONS 

In summary, we have derived the renormalization group equations for a dilute Bose gas. 
To obtain these equations we have to make a distinction between the case of a negative 
and a positive chemical potential. In the latter case we have to take the presence of a con- 
densate into account, and the renormalization group equations are different from the ones 
valid in the unbroken phase. Our philosophy of using the renormalization group approach 
is not to study universal properties of the dilute Bose gas, which are the same as for the 
0(2)— model in three dimensions, but to make quantative predictions about various nonuni- 
versal properties of this system. To achieve this goal one has to find a method to eliminate 
the ultraviolet cutoff dependence inherent to the application of the renormalization group. 
Our knowledge about the two-body scattering problem is sufficient in this respect and we 
can fix the renormalization group equations in this manner. We compared results of the 
renormalization group with mean-field calculations and showed the difference to vanish in 
appropriate limits. We also checked the influence of three-body effects, which turned out 
to be unimportant even in the critical region. Furthermore, we showed that the influence 
of bubble diagrams on the effective interaction can be fairly large, and that the effective 
scattering length a^-'^-^ vanishes when approaching the critical temperature. This confirms 
earlier results PUIP^I and is of importance for a lot of current work concerning condensate 
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properties in trapped Bose gases. It implies that there can in principle be an important 
change in the results of calculations that make use of the nonlinear Schrodinger equation 
when determining these properties. For that purpose the effective scattering length a'^^-^ is 
of importance, and not the two-body scattering length a. 

Next, we derived the renormalization group equations in the broken phase, i.e. for positive 
chemical potential, and showed that the transition to the Bose-Einstein condensed phase 
is contained in this set of renormalization group equations. We found a critical exponent 
u = 0.685, which agrees very well with u = 0.67 found in studies of the 0(2)— model, and 
furthermore calculated the effect of the interactions on the critical temperature of the phase 
transition. The change can be as much as 10 % in the current type of experiments and 
may be measured if one can improve the precision in determining the temperature and the 
central density in the trap. 

As these latter renormalization group equations lead to a dispersion relation which has no 
linear part, i.e. there is no sound mode, we have to restrict ourselves to the regime naA^f^ <^ 1, 
where the linear part of the dispersion is unimportant for determining properties such as 
density and pressure. The region where noA^^ <^ 1 is currently still the most interesting 
one from the experimental point of view although the other region is certainly within reach. 
Including the renormalization of the time derivative and gradient terms in the action is 
expected to resolve the problem of the disappearing of the sound mode, but we postponed 
this to future work. A resolution of this problem would correspond to a resolution of the 
long standing problem of the infrared divergencies in the perturbation expansion around the 
Bogoliubov theory. 

We concluded this work with using the renormalization group to calculate superfluid and 
condensate densities as a function of temperature and also a p — n~^— diagram is presented. 
Of course the renormalization group can be used to find out many more things about the 
dilute Bose gas. Note e.g. that it is in principle also possible to use the renormalization group 
method to find the quantative form of the correlation function as it is the Fourier transform 
of the occupation number A^k- This can also be translated into a differential equation as in 
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Eq. (19) where we now have the distance r as a free parameter. Also the specific heat may 
be calculated. In principle, our results pertain to a homogeneous Bose gas, but in situations 
where the application of a local density approximation is allowed, they are also applicable 
to trapped Bose gases. Moreover, we indicated that it is in principle also possible to set up 
a renormalization group calculation for the inhomogeneous case. 

Finally, we want to note that the procedure of renormalization group as described in 
this article can in principle also be used to study the Kosterlitz-Thouless transition to the 
superfluid phase in a two-dimensional Bose gas. However, there are difficulties in this case 
connected with the fact that all coupling constants are relevant at the critical temperature. 
Nonetheless, work along these lines is in progress since a number of experiments are currently 
under construction which aim at reaching the Kosterlitz-Thouless phase in doubly spin- 
polarized atomic hydrogen adsorbed on a superfluid helium film. 
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In this appendix we go through some of the technicalities of the calculation of the one- 
loop Feynman diagrams for the coupling constants Fn and F12, i.e. essentially the normal 
and anomalous selfenergies, and explicitly show the Hugenholtz-Pines theorem to hold in 
our renormalization group approach. This theorem states that p5| 



where /iSii(O;0) and /;,Si2(0;0) are the irreducible normal and anomalous selfenergies re- 
spectively. In our notation this relation reads 



APPENDIX A 




(A.l) 



Til — ri2 • 



(A.2) 
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Fig. 7 in the main text contains the diagrams renormahzing and in Fig. 12 we depict 
the one-loop diagrams renormahzing Fn. Using the designation of the diagrams as in these 
figures one can show, after going through the combinatorics, that the (infinitesimal) change 
of the vertices after integrating out an infinitesimal momentum shell is given by 

dVn = 4c/F^\ + 4rfFfi + 2dT^[ + AdT^^ + AdVf^ + AdT(^ + ATsd{^) , (A.3a) 



rfFi2 = 2c/Ff2 + 4rfFf2 + AdTg + AdT^^ + 4rfFf2 + 2c/F{'2 + 2F3c/(7^) . (A.3b) 



The last term in both expressions originates from the shift in ao,o required to eliminate the 
linear term dFo(aQo + ao,o) from the action. From Figs. 7 and 12 it is clear that a number 
of diagrams contributing to Fn or F12 are mathematically identical. Most notably we have 
c/Ffi = rfFfi = dVf^ = dTf^: dVf-^ = dVg = dVf^, and dVf^ = dVf^. Therefore, there 
are only six independent diagrams, namely dT^^,dTf^,dT^[,dT^^,dTf^, and dT^2- Using the 
Bogoliubov transformation from Eq. (^), it follows that 

«n«k,n) = ui{blJk,n) + (^-k,-n&-k,-n) , (A.4) 
«na-k,-n) = (ak,na_k,-n) = -Mk^^k( (&k,n V") + (&-k,-n^-k,-n) ) , (A.5) 

with Mk and given in Eq. (p§), and these diagrams are now straightforward to calculate 



by applying the usual Feynman rules ^M. We have 



^^rn = 7:Vodn' = — Vq -K + vl){2N{hujA) + 1) - - e'^'dl, (A.6a) 



1_ . , A^ _ /I , , 1 

2 " 47r"2 " V2^ ■ ' ' 2 



rfFfi = ^Vl{u\vK + ut.vl) {\±^^^^^ + {iN{hujK) \N{huJt.) + 1]) e-^'dl , (A.6b) 



rfrf, = -^Yl {{u\ + 4)^^^^^^ + 1u\vli5N{h^^)\N{hu^) + 1]) e-^Ul, (A.6c) 
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dTf, = -^^Tl + « + vi)PN{hu;j,mhu;j,) + 1]) e~'Ul , (A.6e) 



dT^2 = \yodfi = -^VoUaVa{1 + 2N{hujA))e-''dl , (A.6f) 

introducing dn' = JS?Y.n{ci*K,n^k,n)e~^''dl/2'n'^ and dh = A^J2n{(^l,n(^*-A-n)^~^''dl/2TT'^ with 
which one can easily verify that 

diV^o) = -TT^V = -T^-^r3(2rfn' + dn) . (A.7) 

J- 11 + J- 12 J- 11 + J- 12 

Combining these equations leads to the flow equations for Fn and ri2. However, we here 
concentrate on the Hugenholtz-Pines relation, which is satisfied at / = 0, and therefore 
remains valid if also dTu = dri2 at / = 0. Using the above mentioned equalities of the 
various diagrams, this equality reduces in first instance to 

2dT^^ + dTf^ + T^di,/^) = dVf^ + dVg = dVf^ + dT^^ . (A.8) 

However, due to the f/(l)— symmetry of the action S[ilj*,tp] we have at any value of I that 
r| = uqVq = T12V0 and at / = 0, when Fn = ri2, Eq. ( [A.8|) reduces with the help of Eq. 
O to 



drfi - drfi = Vodn , (A.9) 



or equivalently 



- ^r^« - + = -|^^oUa.a(1 + 2N{nu;A)) . (A.IO) 

Note that, due to the fact that u\ — vj^ = 1, the left-hand side of this equation actually does 
not contain Bogoliubov coherence factors. From the definition of the coherence factors in Eq. 
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(^) it is now easy to show that indeed, at / = 0, we have dTn = dri2 and as a consequence 
that the Hugenholtz-Pines relation remains vahd at any point during the integration. 

Next, we briefly turn to the regime naAff^ <^ 1 where we put = 1 and f k = in 
contributions containing infrared divergencies. In this case, only the first term in the left- 
hand side of Eq. ( [A.IOD contributes, but now with a factor 1 instead of u\. Thus, nothing 
changes with respect to the previous situation, and the Hugenholtz-Pines relation is again 
valid at any point during the renormalization if we still use the exact expressions for u-^ and 

in contributions not containing infrared divergencies, and most importantly in dn. 



APPENDIX B 

In this appendix we discuss the implications of anomalous scaling, focussing on the 
situation at the critical temperature where the coupling constants flow into the fixed point 
(/i*, V^*). We start with the fundamental statement from renormalization group theory that 
the exact n— point vertex function remains identical during renormalization, i.e. in our case 

r(")(p„^i;/i,1/o,T;A) =r(")fe,^,;/x(0,Vb(0,T(/);Ae-') , (B.l) 

where pi < Ae~K The coupling constants obey the derived renormalization group equations 
(without trivial scaling) and T(/) = Te'^K When we are at the critical temperature, we can 
take Ui = and we find after performing the trivial rescaling that 

r(")fe,0;/i,Vo,T,;A) = e-(3-t)'r(")(p,e',0;/i(0,Vb(/),T,(0;A) • (B.2) 

A dimensional analysis then shows that 

r(")(p.e\ 0; M/), Vo{l),m;A) = A'~"^T^-\^, 0; ^, Ml , Ml ^ i) . (b.3) 

Combining these two equations, and taking the limit Z oo in which we approach the fixed 
point, this implies 

rWfe,0;;.,Vo,T.;A)=A3-i(/)rW(-^,0;^,^,oo;l)^ ^'""^^^^'^^^^^^y ) " ^^'^^ 
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The right-hand side has, just hke the left-hand side, to be independent of I, and this leads 
to the conclusion that 

r^'^)*!^) = A5-^(0r(")*(p0 (b.5) 

or that r^"'^*{pi) has to be a homogeneous function of degree (3 — n/2). Thus 

rW*(Ap,) = A=^-trH*(p,) (B.6) 

and we can conclude that anomalous scaling reveals information about the momentum de- 
pendence of the n— point vertex function. In particular, we have e.g. for the selfenergy at 
the critical temperature that 

r(2)*(A;) oc , (B.7) 

and most importantly for the four-point function that 

r(^)*(k, k', K) oc (I k I + I k' I) + a I K I , (B.8) 

which shows that the effective interaction at long wavelengths has to vanish at the critical 
temperature. 

Clearly, the above reasoning is in principle not restricted to the critical temperature and 
a similar argument can be set up also for arbitrary temperatures. As a result, the anomalous 
scaling we find from the set of renormalization group equations derived in Sec. IVA for the 
symmetry broken phase also implies a nontrivial momentum dependence of the coupling 
constants. 
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FIGURE CAPTIONS 



Fig. 1 Typical behavior of the Fourier transform of the interatomic interaction potential 
V"(x-x'). 

Fig. 2 The one-loop Feynman graphs contributing to the renormalization of (a) the chemical 
potential, (b) the two-body interaction, and (c) the three-body interaction. The dot repre- 
sents the two-body vertex Vq and the square the three-body vertex Uq. 

Fig. 3 The p — diagram including (sohd hne) and excluding (dashed line) a three- 
body interaction term for atomic ^^Na at a temperature of O.l/iK. The influence of Uq is 
approximately 1.5 % near the critical density. 

Fig. 4 The ratio of the many-body scattering lengths a^*^ including (solid line) and ex- 
cluding (dashed hne) the bubble diagrams in the flow equations as a function of T/Tc for 
atomic ^^Na at a density of 1.5 10^^ cm~^. 

Fig. 5 Flow diagram resulting from the renormalization group equations. The fixed point is 
indicated with an asterix. 

Fig. 6 The critical degeneracy parameter ricAf^ found from the renormalization group cal- 
culation as a function of a/Afh- 

Fig. 7 The one-loop Feynman diagrams contributing to the renormalization of (a) the linear 
term and (b) the anomalous selfenergy in the broken phase. The filled circle represents the 
vertex Vq and the open circle the vertex T^. 

Fig. 8 The ratio of the effective scattering length a^^^ and the two-body scattering length 
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a for atomic ^^Na at a density of 1.5 10^^ cm~^. The dashed hne corresponds to the result 
from the many-body T-matrix calculation, i.e. a!^^^ — a^^ . The sohd hues correspond to 
the result from the renormalization group approach, i.e. a^^^ = a^*^, with bubble diagrams 
included and excluded. The effect of the bubbles (lower solid hne) is seen to be considerable, 
certainly below the critical temperature. 

Fig. 9 The degeneracy parameter nA^^ as a function of the fugacity C for the ideal gas 
(dotted line), from the many-body T— matrix theory (dashed line) and as found from the 
renormalization group calculation (solid line). 

Fig. 10 (a) The condensate fraction and (b) the superfiuid fraction as a function of tempera- 
ture for a density of 1.5 10^^ sodium atoms per cubic centimeter, both from the many-body 
T— matrix calculation (dashed line) as from the renormalization group calculation (solid 
hne) . 

Fig. 11 The pressure as a function of inverse density for atomic ^^Na at a temperature 
of O.l/iX. The dashed hne is the result from the many-body T— matrix calculation, the 
solid line from the renormalization group calculation. 

Fig. 12 The one-loop Feynman diagrams contributing to the renormalization of the nor- 
mal selfenergy in the broken phase. The filled circle represents the vertex Vq and the open 
circle the vertex F^. 
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Fig. 1 Bijlsma et. al. 




Fig. 3 Bijlsma et. al. 




Fig. 4 Bijlsma et. al. 




Fig. 5 Bijlsma et. al. 




Fig. 6 Bijlsma et. al. 




Fig. 7 Bijlsma et. al. 




Fig. 8 Bijlsma et. al. 
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Fig. 10(b) Bijlsma et. al. 
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